Physics of Spin Casting Dilute Solutions 
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We analyze the spin casting of dilute (ideal) binary mixtures of non-volatile solutes in volatile 
solvents. The system-specific fundamental length and time scales of the process are introduced and 
an analytical description of the thinning of a volatile liquid film simultaneously subject to spinning 
and evaporation is presented for the first time. This reveals spin casting essentially as evaporative 
thinning of a hydrodynamically flattened film. The spatio-temporal composition evolution within 
the thinning film is analyzed. We identify the Sherwood Number as fundamental process parameter 
characterizing whether diffusion (Sh < 1) or evaporation (Sh > 1) dominates the composition 
evolution. Power laws link the process duration, final solute coverage and the spatio-temporal 
evolution of the internal film composition to the control parameters and the properties of the initial 
solution. The analysis is relevant for casting real solutions, in particular for the deposition of 
specifically structured (sub)monolayers ("evaporation-induced self-assembly"). The quantitative 
validity of the analysis can be verified experimentally. Because the analysis provides palpable 
physical insight into the influence of the different process parameters, improvements regarding real 
(non-ideal) conditions can be implemented easily. 

PACS numbers: 68.03.Fg, 47.57.eb, 47.85mb, 82.70.-y 



Introduction. — An excess amount of liquid deposited 
on a spinning, planar, wettable substrate forms a thin- 
ning film of uniform thickness h as the liquid flows out- 
ward (and is spun off) Evaporation of volatile film 
components adds to the film thinning (Fig. |T|). The pro- 
cess of persistent deposition of non-volatile components 
on a spinning support is called "Spin Casting" or "Spin- 
Coating". Coverage, structure, and other properties of 
the final deposit are strongly influenced by the processes 
occurring during film thinning. The nonvolatile com- 
ponents get continuously enriched, eventually exceeding 
saturation. The final deposit evolves from aggregation 
during film thinning in the presence of the solvent. 

The application of spin casting is dominated by de- 
positing polymer films with thicknesses up to tens of 
nanometers Q . But increasingly it is also used to deposit 
structured (sub)monolayers, e.g., specific particle arrays 
( "evaporation- induced self-assembly") [3-ll|. Here we 
study the spin casting of ideal mixtures of a non- volatile 
solute and a volatile solvent to gain universal insights 
into the spatio-temporal composition evolution within 
the thinning film. For the first time the ideal spin cast 
equation is solved analytically. A universal cross-over 
behavior from hydrodynamic- to evaporation-dominated 
film thinning is revealed. The spatio-temporal composi- 
tion evolution within the thinning film is analyzed with 
a diffusion-advection approach. General, palpable rela- 
tions between the process parameters and the composi- 
tional evolution within the film are presented and the 
final solute coverage is calculated. 

The insights are quantitatively valid for dilute solu- 
tions e.g., for a better understanding of evaporation- 
induced structuring of nano deposits. The quantita- 
tive revelation of the spatio-temporal composition evo- 
lution recommends spin casting of diluted solutions for 



early phase 

(spin-off dominated) 



A Z 



flow field '-4b 



spinning support 



angular frequency to 
spin-off coefficient K « to 2 



late phase 

(evaporation dominated) 



solvent evaporation (rate E) 



solute enrichment 



> r 



FIG. 1: Schematics of spin-casting with processes dominating 
early and late stages. 



fundamental studies on nucleation and growth [8|, Il2j . 
Some results, e.g., the predicted solute coverage, are valid 
(semi-)quantitatively even for semi-dilute solutions. In 
addition, the findings improve the understanding of spin 
casting real polymer solutions despite their often indi- 
vidually distinctive, non-ideal rheological and evapora- 
tive behavior, because the physics behind the unveiled 
relations is comprehensive due to the largely analytical 
approach. Therefore our treatment also reveals excellent 
starting points for more detailed, system-specific analyt- 
ical approaches of real systems [16| . 

/. Hydrodynamic- evaporative film thinning - - The 
thinning of a Newtonian, volatile liquid film of thickness 
h on a rotating support is described by (l3j : 



dh/dt= -2Kh 3 -E, 

with spin-off coefficient K = w 2 /(3^), 
speed, v = kinematic viscosity and E = 
The fundamental form of Eq. [TJ 

df/dr = -e-i, 



(1) 

= rotational 
evaporation rate. 
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FIG. 2: Universal thinning curve h/h tr ~ f(t/t sc ) (dashed 
lines: thinning due to spinning/evaporation only). Inset: 
Measured thinning curve in physical units for toluene and 
its fit according to the theory (Eq. [S]). 



is obtained by rescaling £ = h/h tr and r = t/t* c i.e., by 
the system inherent "natural" scales: 



fHr = {E/2K) 1 ' 3 , 
t* sc = (2E 2 K)- 1 /\ 



(3) 
(4) 



h tr is the "transition height" when evaporative and hy- 
drodynamic thinning are equal. t* c is the "reduced pro- 
cess duration" (Eq. 

Eq. [2] is nonlinear and inhomogeneous, but its inverse 
is of zeroth order, and can be integrated: 
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, (5) 



with integration constants so that r = for £ — > oo. Be- 
cause t(£) is bijective, Eq. [5]also yields £(r) respectively 
Mi). 

The duration of the whole process (starting at h — > oo, 
ending at h = 0) is 



t sc = {2n/3 3 / 2 ){2E 2 K)- 1 / 3 = (2tt/3 3 / 2 ) • t* 



(6) 



Fig. [2] shows the universal thinning curve (the inverse 
of Eq. [5]) on scaled axes [ItJ ■ The findings agree with 
experimental results (see inset and [HI)- The transition 
time, ttr, at which h tr is reached, is universal for all spin 
cast processes described by Eq. [1] 



t tr = ((2vr - V31og4) /(4tt)) • i sc a 0.309 • t. 



(7) 



Hydrodynamic film thinning (E — 0, dotted) domi- 
nates «30% of the spin cast time; «70% is evaporation- 
dominated (K = 0, dashed). 

77. Solute concentration evolution — During spin cast- 
ing, non- volatile solute enriches at the free surface (where 
the solvent evaporates) and migrates into the film via 
diffusion. The spatio-temporal evolution of the solute 
concentration c is described by 



d t c = Dd 2 c + Kz 2 (3/i - z) d z c, 



(8) 



with boundary conditions [Ij 



Dd z 



= Ec 



d z c\ z=0 = 0. 



(9) 



The advective term in Eq. [5] is derived from the radial 
velocity field u[r, z) of a Newtonian fluid rotating with 
its solid support (no slip), with a free surface at the top 
(no stress) [l[: 



u(r, z) = 3K r z (h — z/2) . 



(10) 



r and z are the radial respectively vertical coordinates. 
The radial volumetric flux (j) dz is 

4>dz — 2tt r u(r, z) dz — 6n K r 2 z (h — z/2)dz. (11) 

With the continuity equation this yields the thinning- 
induced vertical motion of the horizontal stream lines: 



dZ/dt = -l/(2nr) 



j d r cj>dz' = -K z 2 (3h-z). (12) 
J o 



In order to solve Eq. [8j h(t) is required but not known 
explicitly, only its inverse, t(h) (Eq.EJ. Since h(t) respec- 
tively t(£) are bijective, the time variable can be changed 
to £ d t c = dr/dt ■ d£/dr ■ fl f c = -{2E 2 Ky/ 3 ^ 3 + 1) <%c. 
Rescaling with y = z/h (y 6 [0,1], from substrate to 
surface) avoids the moving boundary. This leads to: 



dec = — 



S7i£ 2 (£ 3 



(^) 2 (3-y) y 

2(e + 1) e 



dyC, (13) 



d v c\ _ = Sh£c\ i; 
y \y—i ^ \y=i 



0, 



(14) 



where the Sherwood number, Sh, parameterizes the ratio 
of evaporative to diffusive mass transport on the charac- 
teristic length scale of the system, h tr - 



Sh = Eh tr /D = E i/3 (2K)- 1/3 /D. 



(15) 



By scaling c with the initial solute concentration cq the 
initial condition is cj^^ = 1. Thus the system is 
parametrized completely by Sh. Eq. [13] can be solved 
numerically with £ as independent variable. Eq. [5] then 
provides t as function of £ i.e., finally, c(t, z). 

III. General aspects of the concentration evolution - 
Based on the solution of Eq. Q2] the impact of Sh (i.e., 
of the individual system parameters K, E, and D) is 
now analyzed. Fig. |3] shows profiles of c during film thin- 
ning for Sh larger and smaller than 1 (i.e., convection 
dominating over diffusion and vice versa) and for film 
heights larger and smaller than h tr , respectively. It re- 
veals the competition between evaporative enrichment, 
spin-off, and diffusive dilution in both, spin-off and evap- 
oration dominated regimes. In general, larger Sh means 
more pronounced gradients in c. If h ^> h tr , solute en- 
richment occurs only locally near the free surface and 
c ~ cq near the substrate. If h <C h tr , c also increases 
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FIG. 3: Solute concentration profiles (resulting from Eq. I13p 
for three different Sherwood numbers 9i at three different 
moments respectively heights h/h tr - 



near the substrate. Sh = 1 and h = h tr mark the transi- 
tion: The solute gradient just reaches the film/substrate 
interface and c just begins to increase globally. 

Fig. |4] shows as function of h/htr- A) The total solute 
amount N per unit area A (N/A = f Q cdz); B) c/co 
at the surface respectively substrate/film interface; C) 
The difference between the surface and the substrate / film 
interface concentrations, Ac/cq i.e., the relative enrich- 
ment. Both, h = htr and Sh = 1 mark transitions be- 
tween distinctly different behaviors. I) For h 3> htr spin- 
off dominates film thinning. Therefore, globally c/co re- 
mains approximately constant and N/A decreases. Nev- 
ertheless, there is a surficial solute enrichment due to 
evaporation. But this Ac/co only becomes substantial 
for Sh > 1. It has a maximum at h = h pea k (h pea k w h tr 
for Sh = 0(1)). II) For h <C h tri evaporation domi- 
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FIG. 4: For various Sherwood numbers Sh are shown as func- 
tion of the reduced film thickness h/htr- A) Total solute 
amount N per per unit area A (scaled by cohtr)', B) Con- 
centrations c at the surface respectively film/substrate inter- 
face (scaled by Co); C) Ac, the difference between the surface 
and the substrate/film interface concentrations, respectively 
(scaled by Co). 
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FIG. 5: Normalized final solute coverage, N(h — > 0)/(Acohtr) 
and c/co at the film surface for h = htr as function Sh. 

nates. Solvent loss leads to increasing c/cq while N/A 
and Ac/ cq remain constant (because spin-off becomes 
negligible). Remarkably, for h < h tr , N/(Acoh tr ) re- 
mains approximately constant and independent from Sh 
whereas Ac/cq increases with increasing Sh. 

Fig. [S] shows c(z — h — h tr )/co, the concentration at 
the free surface for h — h tr , and N(h — > 0)/(Acoh tr ), 
the rescaled final coverage. Both are plotted as func- 
tion of Sh. For Sh <C 1 diffusion dominates and c is 
mostly homogeneous (Fig. [3]). In this case coEt sc is the 
total amount of solute that is not spun off with the 
solvent. Distributing this into a film with ht r yields 
c/c = Et sc /htr = 2tt/3 3 / 2 fa 1.2 (Fig. \5§. The final 
solute coverage for Sh < 1 is (in agreement with [3]) 

T = N(h -> 0)/A w c h tr ~ 0.8 c {K/E)- 1 ' 3 . (16) 

For Sh ^> 1 it is only w 6% lower. 

IK Relative surficial enrichment maximum — Fig. [6] 
shows Acmax/co and its spatio-temporal properties as 
function of Sh. Symbols denote the results from the nu- 
merical analysis. The solid lines show power laws for 




FIG. 6: Acmax/co as a function of Sh and related spatio- 
temporal properties (h vea k=h\ Ac and f P eafe=i| Ac )• 
The symbols indicate numerical results, the lines the under- 
lying power laws for Sh < 1. 
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Sh < I which are rationalized by analyzing the underly- 
ing processes. 

Panel A) shows h peak rescaled by h tr i.e., the film thick- 
ness at which Ac = Ac max as function of Sh. Ac max 
emerges from the competition between evaporative en- 
richment, spin-off, and diffusional equilibration. Evap- 
oration and spin-off dominate at opposite ranges of h. 
For large ft, surficial spin-off efficiently suppresses en- 
richment: Ac,nax requires 2Kh peak < E (from Eq. Q]). 
However, for small ft, diffusional equilibration is fast, so 
Ac m ax requires D/h pea k < E. Optimizing both con- 
ditions simultaneously {E is linear and protagonistic in 
both cases, and therefore cancels out) reveals the same 
power law as numerics (which supplies the prefactor): 

h P eak « {D/K) 1 / i ps 1.2 htr Sh' 1 '* . (17) 

Panel B) shows the time t pea k rescaled by t sc i.e., the 
time at which Ac = Ac max . Before reaching h pea k , thin- 
ning is dominated by spin-off. Hence t pea k can be esti- 
mated by inserting h pea k into h = (2Kt)~ 1 / 2 (the solution 
to Eq.CQ with E = 0): 

t P eak/tsc « 0.31 Ef/'K-V'D-V' w 0.35 VSh. (18) 

Panel C presents Ac max /co, reflecting the balance be- 
tween evaporative enrichment and diffusive equilibration 
i.e., Sh\ h ^ = Eh peak /D = Sh 3/ \ With Eq. El this 
means: 

Ac ma x/c w 0.46 E K~ 1/4 D- 3/i ps 0.55 Sh 3/ \ (19) 

Discussion and Conclusions — Section I introduces 
the system-specific fundamental length and time scales 
(h tr , t* sc ) for the spin casting process of an ideal Newto- 
nian volatile liquid. These reduce the general spin cast 
equation (Eq.QJ to its fundamental form. The equation is 
solved analytically for the first time revealing a universal 
film thinning behavior. The total spin coating time, t sc is 
calculated (Eq. [BJ as function of the process parameters 
(E, K, D). For any combination of these parameters, in 
the first 30% of the process time, thinning is governed 
by hydrodynamics. During the residual 70% of the time 
evaporative thinning dominates until complete drying. 

Sections II through IV analyze the spin cast process 
assuming an ideal mixture of a volatile solvent and a 
non-volatile solute. It is found that the spatio-temporal 
evolution of the solute concentration within the thinning 
film is universally characterized by its Sherwood number, 
Sh, which reflects the competition between evaporative 
solute enrichment and diffusional dilution on the system- 
inherent fundamental length scale, htr- For Sh < 1 (dif- 
fusion dominates) the spatio-temporal occurrences of the 
solute concentration maxima are related to the process 
parameters via universal power laws (Eqs.EIHHJ andHH 
see Fig. |5|) . These findings are rationalized semiquanti- 
tatively with the underlying physics. At last, the final 



solute coverage Y is calculated from the process parame- 
ters (Eq.QU). 

To examine the relevance of our analysis for real cases, 
where E and K are not constant and depend on the solute 
concentration (E = E(c), K = K(uj, z/)), we assume a 
typical solvent (e.g. toluene) with molar mass M s ps 
O.lkg/mol, density p s ss 10 3 kg/m 3 , v — 6 ■ 10~ 7 m _2 s _1 , 
and E = 10~ 6 m/s [18]. We consider two examples: A) A 
typical polymer solution 0; B) A nanoparticle solution 
for submonolayer deposition. 

Case A): Polymer with M P = lOOkg/mol, p P = 
10 3 kgm~ 3 , and co, P = lOkgm -3 (i.e., mole fractions 
x , P = 10~ 6 and x , s ps 1), D P = 10" 11 n^s -1 and 
K = 10 11 m _2 s _1 (~ 4000rpm). This yields: h tr ~ 
1.7-10- 6 m (Eq.E), Sh w 0.17 (Eq.HH), c P (h tr )/ca^ « 1.2 
(Figs. [3] and EJ and T p w 17 * lO^kgm- 2 (Eq. [T6J) i.e., 
a final film thickness of ps 17 nm. This is reasonable Q. 

Case B): Spheres with radius r NP — 36 nm. _D NP = 
10 -11 m 2 s _1 (estimation: Stokes-Einstein) as in case A) 
and thus ht r and Sh are identical to case A). A 50% 
monolayer coverage means r NP ps 10 m -2 i.e., co, NP ~ 
5 • 10- 5 mol/m 3 (x , NP ~ 5 • 10" 9 , x 0:S ps 1) with Eq.HH 

In both cases cq is small (xo,s *C l,xo, NP <C 1) i.e., 
the initial solution behaves approximately ideally. This 
means film thinning to ht r will follow Eq. [T] with K and 
E remaining approximately constant because thinning is 
hydrodynamically dominated and the solute concentra- 
tion barely changes. Hence Y can be calculated with 
Eq. [16] because it is not affected by changes in K and E 
occurring after reaching /i tr . 

For both examples Sh < 1. Therefore the data de- 
picted for Sh < 1 in Figures [5] and [6] are relevant because 
h P eak, tpeak, and Ac max occur at h > h pea k when the so- 
lutions still behave approximately ideally. Accordingly, 
also Eqs. El [18] and [19] are applicable. 

Of course, evaporative film thinning at ft < h pea k 
will eventually increase the solute concentration and 
thus change E(c) and K(uj,v). However, in particular 
for "evaporation-structured" submonolayer particle ar- 
ray deposition, non-ideality will only become relevant for 
h <C h tr , when h is already in the range of r NP , because 
typically a; , NP < 1- 

Our analysis provides valuable starting-points for im- 
proved approaches (other aspects) beyond the idealized 
continuum hydrodynamic-evaporative model. For in- 
stance, for low solubilities our analysis reveals where so- 
lute aggregation occurs first : A) at the top surface due 
to sufficiently high evaporative solute enrichment or B) 
globally homogeneous (or heterogeneously at the sub- 
strate) because diffusion keeps the (increasing) concen- 
tration approximately uniform. 

Last not least, our analysis predicts a quantifiable be- 
havior based on measurable parameters of the initial so- 
lution. In particular the predicted Y and h(t) [18] are 
easily measurable. Thus not only the validity of the ap- 
proach can be evaluated. It is possible to specifically ad- 
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dress the influence of individual parameters of the multi- 
parameter spin casting process. For instance, parame- 
ters not related to mixing properties (e.g., temperature 
changes induced by evaporation) can be probed, because 
a sufficient reduction of cq will per definition render the 
solution behaving approximately "ideally" . All this pro- 
poses our approach as an excellent basis for incorporating 
specific (non)-linear properties of specific systems. Due 
to the predominantly analytic approach the basis still 
will be of general character and reveal palpable physical 
insights. 
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